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Abstract: Comprehensive assessments of ecosystem services in environments under the influences of 
human activities and climate change are critical for sustainable regional ecosystem management. Therefore, 
integrated interdisciplinary modelling has become a major focus of ecosystem service assessment. In this 
study, we established a model that integrates land use/cover change (LUCC), climate change, and water 
retention services to evaluate the spatial and temporal variations of water retention services in the Loess 
Plateau of China in the historical period (2000—2015) and in the future (2020—2050). An improved 
Markov-Cellular Automata (Markov-CA) model was used to simulate land use/land cover patterns, and 
ArcGIS 10.2 software was used to simulate and assess water retention services from 2000 to 2050 under 
six combined scenarios, including three land use/land cover scenarios (historical scenario (HS), ecological 
protection scenario (EPS), and urban expansion scenario (UES)) and two climate change scenarios 
(RCP4.5 and RCP8.5, where RCP is the representative concentration pathway). LUCCs in the historical 
period (2000-2015) and in the future (2020-2050) are dominated by transformations among agricultural 
land, urban land and grassland. Urban land under UES increased significantly by 0.63X103 km?/a, which 
was higher than the increase of urban land under HS and EPS. In the Loess Plateau, water yield decreased 
by 17.20X10° mm and water retention increased by 0.0910° mm in the historical period (2000-2015), 
especially in the Interior drainage zone and its surrounding areas. In the future (2020-2050), the pixel 
means of water yield is higher under RCP4.5 scenario (96.63 mm) than under RCP8.5 scenario (95.46 
mm), and the pixel means of water retention is higher under RCP4.5 scenario (1.95 mm) than under 
RCP8.5 scenario (1.38 mm). RCP4.5-EPS shows the highest total water retention capacity on the plateau 
scale among the six combined scenarios, with the value of 1.27X10° mm. Ecological restoration projects 
in the Loess Plateau have enhanced soil and water retention. However, more attention needs to be paid 
not only to the simultaneous increase in water retention services and evapotranspiration but also to the 
type and layout of restored vegetation. Furthermore, urbanization needs to be controlled to prevent 
uncontrollable LUCCs and climate change. Our findings provide reference data for the regional water and 
land resources management and the sustainable development of socio-ecological systems in the Loess 
Plateau under LUCC and climate change scenarios. 
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1 Introduction 


Ecosystem services refer to the benefits that humans derive directly or indirectly from ecosystems 
(Costanza et al., 1997). Ecosystem services assessment frameworks such as Millennium 
Ecosystem Assessment and Economics of Ecosystems and Biodiversity (TEEB) Programs 
emphasizing the close relationship between ecosystems and human well-being, have greatly 
advanced ecosystem services research (Perrings et al., 2011). With the application of ecosystem 
services research on a global scale, numerous countries and government organizations have 
utilized it in the formulation of development policies and the rational allocation methods of 
natural resources, which are projected to be the long-term trends in ecosystem services research 
pertaining to the global sustainable development (Bai et al., 2019). 

Water retention services are defined as the capacity of an ecosystem to intercept and store 
precipitation, which reflects the sustainability of regional water resources to a certain extent (Bai 
et al., 2011; Quintas-Soriano et al., 2014). Water retention services are key factors determining 
the continuous supply of regional water resources and regional development (Feng et al., 2016a; 
Hackbart et al., 2017). On a large spatial-temporal scale, land use and climate change are the 
major factors inducing the changes in water retention services (Su and Fu, 2013; van Vliet, 2019), 
however, the processes and methods of these influences are complex, especially in the future, 
with uncertainties. A large number of studies have attempted to investigate the mechanisms 
underlying these relationships and quantify the influencing processes, in order to support regional 
ecological construction and sustainable development (He et al., 2019; Wen and Théau, 2020; Hu 
et al., 2021; Wang et al., 2021). 

Since the 20" century, large-scale urban expansion and vegetation restoration have been the 
main land use/cover changes (LUCCs) in China. Over the past 30 a, the Loess Plateau in China 
underwent vigorous urbanization, with an urbanization rate of approximately 30% in the whole 
Loess Plateau, while the urbanization rate of cities is about 80% such as Xi'an (Wu et al., 2019). 
Studies have shown that rapid urbanization, a global phenomenon, has significantly changed the 
relationship among atmosphere, hydrosphere, and biosphere, leading to a decrease in regional 
water storage and the expansion of rural settlements (Li et al., 2016). Furthermore, the expansion 
of construction land that performs at the expense of ecological lands, such as cultivated land, forest 
land, and grassland (Peng et al., 2016), is the direct cause of the decline in ecosystem services such 
as water yield. However, population accumulation and economic development are the indirect 
causes of the decline in ecosystem services (Peng et al., 2017). Meanwhile, a series of vegetation 
restoration and protection measures have been taken for promoting soil-water retention services in 
the Loess Plateau, such as the Grain for Green Project (GGP), the Natural Forest Protection Project 
(NFPP), and the Three-North Shelter Forest Program (TNSFP), which have markedly increased the 
vegetation cover in the region (Chen et al., 2015). Su and Fu (2013) pointed out that during 
1975-2008 in the Loess Plateau, the areas of woodland and grassland increased by 0.47x10* and 
1.90x10* km’, respectively, while the area of farmland decreased by 2.5x10* km’. Especially after 
the implementation of the GGP (2000-2008), a large area of steep slope farmland was converted 
into grassland and woodland. On the other hand, vegetation restoration exerts negative effects on 
water retention as well. For instance, LUCC caused by the GGP has decreased water yield by 38% 
in the Loess Plateau (Feng et al., 2012). Furthermore, 0.09x10° km? land that changed following 
the implementation of the GGP during 2000-2015, has resulted in a decrease of runoff by 1.5x10? 
m?in the Loess Plateau (Deng et al., 2019). Therefore, vegetation restoration requires appropriate 
trade-offs and optimized management in the Loess Plateau. 

Water retention services are directly affected by precipitation in the Loess Plateau and tend to 
be more sensitive to climate change than those in other global regions (Jiang et al., 2018). 
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Numerous studies have attempted to reveal the significant changes of water retention services in 
the Loess Plateau, a region that has been under a warming and drying climate trend over the past 
30 a (e.g., Zhou et al., 2005; Su and Fu, 2013; Fu et al., 2017). Feng et al. (2016b) pointed out 
that the average runoff decreased gradually from 1961 to 2009 at a rate of 0.9 mm/a in the Loess 
Plateau; they attributed this phenomenon to a decrease of precipitation. Fu et al. (2017) reported 
that the decrease of precipitation is the key factor influencing the decrease of the Yellow River 
runoff in the Loess Plateau. The Loess Plateau occupies a vast area, so its key elements such as 
land use and climate are heterogeneous in geographic space. Clarifying the spatial characteristics 
and trends of land use and climate change and quantifying their impacts on water retention 
services are indispensable for reviewing and summarizing the suitability of existing development 
models and exploring future sustainable development paths in the Loess Plateau. 

Spatially explicit modelling under scenarios is a useful method for quantifying the impacts of 
land use and climate change on ecosystem services (Martnez-Harms and Balvanera, 2012; Kim et 
al., 2013; Hu et al., 2020). Among modelling tools of ecosystem services, the integrated valuation 
of ecosystem services and tradeoffs (INVEST) model is the most commonly used approach for 
reflecting the impacts of land use and climate change on hydrological processes at multi-scales. 
However, the model is used as a "black box" due to non-spatial input parameters, which simplify 
the evaluation process to a certain extent, yet bringing difficulty to parameter verification (Ochoa 
and Urbina-Cardona, 2017; Chen et al., 2018; Kim and Jung, 2020). Spatial input parameters 
could provide better feedback than non-spatial input parameters, considering the spatial 
variability of climatic and geographic attributes at large scales in the Loess Plateau. Spatially 
explicit modelling using Markov-Cellular Automata (Markov-CA) model is increasingly being 
applied for multi-scenario spatial-temporal simulations of future land use at regional scales (Du et 
al., 2020). In this study, the future land use demand was calculated using Markov model, which is 
based on the initial state probability matrix formed by land use in two periods (Muller and 
Middleton, 1994; Guan et al., 2011). Future land use simulation (FLUS) model, a typical CA 
model, was used for reflecting land supply and differences in actual land use changes. FLUS 
model integrates artificial neural network (ANN) algorithm and roulette wheel selection (RWS) 
mechanism (Liu et al., 2017) and is extensively used in multi-scenario land use simulations with 
high accuracy (Liang et al., 2018; Lin et al., 2020). 

In this study, we aimed to comprehensively understand the nexus relationship between land use 
and water retention services in the Loess Plateau using integrated modelling methods. Markov 
model was first modified to simulate multi-scenario land use demands and then integrated with 
ANN-based suitability probability estimation (ANN-SPE) and self-adaptive inertia and 
competition mechanism (SAICM) modules of FLUS model to calculate the probabilities of 
occurrence for different land use types and their pattern variations. Subsequently, a spatially 
explicit water retention services model was developed using ArcGIS 10.2 (ESRI, Redlands, 
California, USA). 

More specifically, the main objectives of this study are to: (1) develop an integrated modelling 
framework, including modified Markov-FLUS, water retention assessment model, and 
multi-source datasets via various mapping techniques; (2) simulate the spatial-temporal variations 
in land use types and water retention services in the Loess Plateau in the historical period 
(2000-2015); (3) define land use scenarios of historical scenario (HS), ecological protection 
scenario (EPS), and urban expansion scenario (UES), and obtain climate products under RCP4.5 
and RCP8.5 (RCP, representative concentration pathway)scenarios (IPCC, 2017) in the simulation 
period (2020-2050); (4) determine the spatial-temporal variations in water retention services 
based on different combinations of aforementioned land use and RCP scenarios; and (5) 
quantitatively examine the complex influence of land use and climate change on water retention 
services in the Loess Plateau. This study furthers our understanding of the relationship between 
local water and land resources and the results can provide scientific support for regional 
ecosystem management in the Loess Plateau. 


202205.00021v1 


chinaXiv 


ChinaXivA F #AF! 


SUN Dingzhao et al.: Scenario simulation of water retention services under... 


2 Materials and methods 


2.1 Study area 


The Loess Plateau (33°41'—41°16'N and 100°52'—114°33'E) is the largest (approximately 
0.64x10° km?) and most concentrated typical loess landform in the world (Fig. 1). It is 
characterized by a temperate continental monsoon climate, with annual average temperature and 
annual reference evapotranspiration (ETo) of 3.6°C—14.3°C and 1060.00 mm, respectively (Li et 
al., 2012). The region experiences little and highly spatially heterogeneous precipitation that 
mainly occurs from June to September during rainstorms (Zhang et al., 2014), accounting for 
60%-70% of the total annual precipitation in the Loess Plateau and contributing 90% of the 
annual total sediment load of the Yellow River (Fu et al., 2017). 
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Fig. 1 Overview of the Loess Plateau 


In the Loess Plateau, increasing levels of human disturbance (population increasing from less 
than 40x10° in 1949 to over 100x10° in 2000; Fu et al., 2017) and continuous climate change 
have led to serious ecological problems. For example, per capita water availability of the Loess 
Plateau is only one-fifth of the national average in China. Over-urbanization, over-cultivation, 
agricultural land abandonment, and vegetation destruction have also accelerated ecological 
degradation of the Loess Plateau, making it an ecologically fragile area with the most severe 
water and soil loss in the world (Wang et al., 2018). Since the 1980s, a series of ecological 
restoration projects (e.g., the TNSFP, NFPP, and GGP) have been initiated in the Loess Plateau, 
with the aim of alleviating water and soil loss. In 1999-2017, vegetation cover increased by 
nearly 33% (Lü et al., 2012), however, the negative effects of increased vegetation cover are 
creating new challenges for regional water resources management in the Loess Plateau. Several 
studies have pointed out that water demand for increased vegetation cover and canopy ETo caused 
by vegetation restoration measures are putting pressure on limited water resources in the Loess 
Plateau (Sun et al., 2006; Feng et al., 2016a; Shao et al., 2019; Wang et al., 2021). In addition, the 
large-scale water consumption associated with large-scale urban expansion cannot be ignored, and 
water resources management activities have become more complex in the context of general 
temperature rise and precipitation decline in the Loess Plateau (Su and Fu, 2013). Consequently, 
long-term explicit integrated assessments of water resources in the Loess Plateau are required. 
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2.2 Methods 


Modelling was carried out using the following four steps in this study (Fig. 2). First, summarizing 
scenarios (HS, EPS, and UES) and potential driving factors (natural environment, socio-economic 
status, and relevant planning policies in the Loess Plateau) of LUCC as background knowledge. 
Second, adopting an improved Markov model to predict future land use demands under multiple 
scenarios (HS, EPS, and UES) based on land use status in the historical period (2000-2015) of the 
Loess Plateau; utilizing a logistic regression model to evaluate and select factors that potentially 
drive LUCC; and inputting the driving factors and land use demands into FLUS model to simulate 
the potential spatial patterns of land use in the future from 2020 to 2050. Third, calculating actual 
evapotranspiration by integrating soil physical-chemical properties, meteorological data, 
topographic indices, leaf area index (LAI), and Zhang coefficient (Donohue et al., 2012); and 
deriving precipitation datasets according to RCP scenarios. Fourth, utilizing Budyko hypothesis to 
integrate the abovementioned data in the second and third steps to obtain the long-term evaluation 
results of water retention services in the historical period (2000-2015) and the multi-scenario 
simulation results of water retention services in the future (2020—2050) in the Loess Plateau. 


Natural environment Future planning policy 


Land use in the 
historical period 


LUCC potential 
impact factors 


Zhang Leaf Topographic 
coefficient area index index 


Soil eee ee ae | 
eee ee ae | 


Water retention service dynamic in the Water retention service dynamic in the 
future (2020-2050 


Fig. 2 Integrated modelling framework in this study. HS, historical scenario; EPS, ecological protection 
scenario; UES, urban expansion scenario; LUCC, land use/cover change; FLUS, future land use simulation; RCP, 
representative concentration pathway. 


2.2.1 Data sources 


Multi-source data used to assess water retention services in the Loess Plateau are shown in Table 
S1. Among them, land use/cover dataset was processed using the European Space Agency (ESA) 
(https://www.esa-landcover-cci.org) product classification scheme. The ESA climate change 
initiative dataset, which uses machine learning and unsupervised algorithms, was able to fuse 
3.0x10° consecutive land cover maps, and the accuracy of the dataset has been verified on a global 
scale (Li et al., 2018). We reclassified the land use/cover dataset into 10 categories (agricultural 
land, evergreen broad-leaved forest, deciduous broad-leaved forest, evergreen coniferous forest, 
closed forest, grassland, shrubland, urban land, water bodies, and others) according to the principles 
of land iteration efficiency and the data used to support local land management planning; in 
particular, the area of deciduous coniferous forest was small, therefore, it was classified as closed 
forest. Digital elevation model (DEM) and the spatial distributions of population and GDP (gross 
domestic product) were derived from the Resource and Environment Data Cloud Platform, the 
Chinese Academy of Sciences (http://www.resdc.cn/Default.aspx). Soil attribute data were from the 
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Environmental and Ecological Science Data Center for West China (http://westdc.westgis.ac.cn/). 
Precipitation data (2000-2015) and other meteorological data (including temperature, relative 
humidity, barometric pressure, wind direction speed, and sunshine hours) during 1950-2015, which 
were mainly used to calculate ETo, were obtained from the China Meteorological Data Service 
Center (http://data.cma.cn/), and their spatial distributions were derived using ordinary kriging 
interpolation method. LAI dataset was extracted from Land-Atmosphere Interaction Research 
Group at Sun Yat-sen University, China (http://globalchange.bnu.edu.cn/). Future climate data were 
generated on the basis of the National Aeronautics and Space Administration (NASA) Earth 
Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset. NEX-GDDP couples the 
measured values of global sites and data of 21 types of general circulation models (GCMs) in 
Coupled Model Inter Comparison Project Phase 5 (CMIP5); the dataset includes three climatic 
variables: daily precipitation, daily maximum temperature, and daily minimum temperature, for the 
periods of 1950-2005 (historical run) and 2006-2100 (RCP4.5 and RCP8.5 runs) (van Vuuren et al., 
2011). Studies have demonstrated that NEX-GDDP can successfully reproduce the spatial patterns 
of extreme precipitation in China (Chen et al., 2017). Specifically, RCP4.5 and RCP8.5 represent an 
intermediate scenario and an extreme scenario where policy interventions maintain greenhouse gas 
emissions below 4.5 and 8.5 W/m‘, respectively. Under RCP4.5 and RCP8.5 scenarios, the global 
average surface temperature is predicted to increase by 0.9°C—2.0°C and by 1.4°C—2.6°C, 
respectively (Meinshausen et al., 2011). The grid data were downscaled to a global spatial resolution 
of 0.25°x0.25° (Raghavan et al, 2018; Sahany et al, 2019) and projected under 
"Krasovsky_1940_Albers"” with a resampled resolution of 1 km. All the data mentioned above were 
mainly processed using ArcGIS 10.2 software. 

2.2.2 Scenario designs 

Conflicts between agricultural land protection and urbanization are major challenges to land use 
management in the Loess Plateau. Rapid urbanization promotes intensive socio-economic 
development; however, large areas of agricultural land and unused land are occupied, leading to 
potential challenges such as rural hollowing (Fu et al., 2018). Additionally, ecological restoration 
projects may reduce available water resources (Liang et al., 2015). 

To clarify these complex trade-off relationships mentioned above, we designed three typical 
LUCC scenarios based on the interactions among agricultural land, urban land, closed forest, and 
grassland under the influence of climate change in the future (2020-2050) in the Loess Plateau. 
Scenario 1 (HS) is that the variation trends and rates of different land use types are the same as 
those in the historical period (2000-2015). Scenario 2 (EPS) moderately slows the degradation of 
agricultural land, expediates the expansion of forest land and grassland under the closure 
management, and limits the expansion of urban land; specially, the coverage of closed forest and 
grassland increases 1.2 times faster while agricultural land decreases 0.8 times more slowly than 
those under HS. Scenario 3 (UES) is that urban land expansion is rapid and becomes uncontrollable, 
resulting in less effective ecological restoration under this scenario; specially, the rates of urban land 
expansion and vegetation coverage increase are 1.2 and 0.9 times of those under HS, respectively. 

Because of the difficulty in reducing greenhouse gas emissions with the current socio-economic 
development condition (Okkan and Kirdemir, 2016), RCP4.5 and RCP8.5 were used to simulate the 
potential effects of LUCC on water retention services under current and extreme greenhouse gas 
emission scenarios in the Loess Plateau. 

Overall, we created six scenarios by combining land use scenarios (HS, EPS, and UES) and 
climate change scenarios (RCP4.5 and RCP8.5): RCP4.5-HS, RCP4.5-EPS, RCP4.5-UES, 
RCP8.5-HS, RCP8.5-EPS, and RCP8.5-UES. 

2.2.3 LUCC simulation 

(1) Improved Markov-CA model. The prediction results of Markov-CA model are only related to the 
conversion probability matrix of two periods of land use maps, which have no aftereffect (Guan et 
al., 2011). Scenario weight matrix based on the initial state transition probability which reflects the 
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conversion intensities of different land use types relative to their initial states in specific policy 
scenarios, was introduced to control the multi-scenario expansion simulation capabilities of 
different land use types. 

The improved Markov-CA equations are as follows: 
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where Pin) is the transition probability matrix of land use types; n is the number of land use types; P, 
is the modified matrix of the initial state probability matrix Pj after the introduction of scenario 
weight matrix to change the weights of different land use types; and P; is the standardized P,,, 
which ensures the sum of the conversion probabilities of Equation 3. 

(2) FLUS model. We entered the selected land use driving factors using logistic regression and 
ROC (receiver operating characteristic; its value indicates a better effect when it is greater than 0.7) 
analysis (Eq. 4; Pontius and Schneider, 2001) into ANN-SPE model to calculate the probabilities of 
occurrence for different land use types. Afterward, we adopted SAICM model to improve the 
randomness of LUCC and increase the accuracy of land use pattern simulations (Liu et al., 2017; 
Liang et al., 2018). The parameterization of FLUS in the Loess Plateau has been completed in the 
existing study (Sun and Liang, 2021). Then, we used overall accuracy (OA) and Kappa coefficient 
to assess the accuracy of land use pattern simulations (Pontius and Schneider, 2001; Rwanga and 
Ndambuki, 2017), where the situation of Kappa coefficient value higher than 0.8 and OA value 
higher than 0.8 indicates a reasonably high accuracy (Pontius, 2000; Sun and Liang, 2021). LUCC 
simulations were completed when land use demands are met. 
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where P; denotes the probability that land use type i is present in the grid of interest; æ is the 
intercept; 2; is the regression coefficient of the i® driving factor; and x; is the i® driving factor. 

(3) Trade-off. We utilized partial correlation analysis to quantitatively explain the trade-off or 
synergic relationships of specific LUCCs by SPSS software (International Business Machines 
Corporation, Armonk, New York, United States) (Eqs. 5 and 6). The correlation coefficient values 
of higher than, lower than, and equalling to zero indicate that variables exhibit trade-off 
relationship, synergic relationship, and no relationship, respectively (Howe et al., 2014). All the 
results were subjected to the t-tests. 
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where rjis the correlation coefficient between land use type i and land use type j; y is the year; 
Lya and Lyo) represent the area of land use type i and land use type j in a particular year y (km°), 


respectively; L and L; represent the average area of land use type i and land use type j in the 


study period (km?), respectively; rjx is the partial correlation coefficient between land use type i 
and land use type j when the influence of land use type k is controlled; rix is the correlation 
coefficient between land use type i and land use type k; and rj is the correlation coefficient 
between land use type j and land use type k . The eighth-order partial correlation coefficient 
between land use type i and land use type j was progressively obtained when the other eight land 
use types were simultaneously controlled. 

2.2.4 Water retention services assessment 

Based on Budyko hypothesis (Donohue et al., 2012) and annual precipitation data, we calculated 
regional water yield (mm), i.e., the difference between precipitation and AET (transpiration by 
vegetation and evaporation from the ground surface) (Sun et al., 2006), as follows: 

Yield = P—AET , (7) 
where P is the annual precipitation (mm); and AET is the actual evapotranspiration (mm), which 
was calculated from the original parameters such as ETo, LAI, and Zhang coefficient. Especially, the 
annual precipitation amount and frequency (N) data from 133 weather stations in the study area 
were collected, the annual average Z-value for each station was calculated (Z=0.2N; Donohue et al., 
2012), and the ordinary kriging interpolation was performed. 

Water retention (mm) was calculated using water yield, topography, soil retention properties, and 
soil texture, as follows (Wang et al., 2019): 


Water retention = min (1 2) xmin (1 eee =) xmin| 1, ust x Yield, (8) 
V 3 300 


where V is the dimensionless flow velocity coefficient that reflects the water flow resistance; TI is 
the topographic index; and Ksatis the soil saturated hydraulic conductivity (cm/d) calculated with the 
help of Neuro Theta (Minasny and McBratney, 2002) using sand, silt, and clay contents of soil. 


3 Results 


3.1 Land use simulations and trade-off analysis 


3.1.1 Spatial-temporal changes of land use types 


ROC analysis revealed that the driving factors (slope, aspect, precipitation, temperature, 
population, GDP, distance to cities, roads, railways, and rivers) (Sun and Liang, 2021) in all 
logistic equations greatly explained the variations in LUCCs, and that evergreen broad-leaved 
forest and agricultural land had the highest (0.962) and lowest (0.738) ROC values, respectively. 
The land use maps of 2011—2015 were used to analyze the simulation accuracy of FLUS model, 
and Kappa coefficient and OA values were in the ranges of 0.801—0.821 and 0.813-0.886, 
respectively. These results suggested that LUCC simulation with the calibrated LUCC model was 
spatially consistent and accurate. 

The overall patterns in the major LUCCs are relatively consistent under different scenarios in 
the future (2020-2050), however, distinct variations are observed in several regions. Under UES, 
the rapid expansion of urban land (0.63x10° km?/a) mainly occurs in the southern part, and the 
reduction of agricultural land (0.58x10? km?/a) is also noticed. Under EPS, grassland increases 
significantly (0.20x10? km?/a), while deciduous c, evergreen coniferous forest, and closed forest 
all expand to varying degrees. The LUCC under HS shows a medium situation between the two 
scenarios mentioned above. Overall, agricultural land is predicted to shrink by 5.55%—7.39%, and 
is mainly observed in the central hilly-gully regions (Fig. 3), while urban land is predicted to 
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increase by 1.34—2.20 times. However, the complex feedback relationship between urbanization 
and socio-ecological system still requires further investigation in the Loess Plateau. 


(a) HS in 2020 N (b)HSin 2030 


4P 


(c)HSin2040 
G 


(©) EPS in 2020 _ 


s 


(D EPS in 2030 


G 


(g) EPS in 2040 _ 


Cs 


(h) EPS in 2050 _ 


¥ 


(i) UES in 2020 


K 


(j) UES in 2030 


A 


050 _ 
aoa 


Legend 
Land use type 


(Agricultural land mm Deciduous broad-leaved forest EE Closed forest EE) Shrubland mmi Water bodies 
mm Evergreen broad-leaved forest Evergreen coniferous forest ME Grassland mE Urban land æm Others 


Fig. 3. Simulations of land use type patterns under historical scenario (HS), ecological protection scenario (EPS), 
and urban expansion scenario (UES) in the Loess Plateau in 2020 (a, e, i), 2030 (b, f, j), 2040 (c, g, k), and 2050 
(d, h, 1) 


3.1.2 Land use conversion and trade-offs 


The spatial-temporal variations of land use types shown in Figure 4 demonstrate frequent 
conversions among agricultural land, urban land, and grassland in the future (2020-2050). 
Conversions from agricultural land to urban land in the south of the Loess Plateau mainly occur in 
existing cities and provincial capitals, especially under UES. Areas of agricultural land growth 
intersect with areas of agricultural land reduction under HS. However, under EPS, local and 
concentrated growth or shrinkage of agricultural land is noted. Under UES, significant reductions 
of agricultural land are observed locally, especially the conversion to urban land. Under EPS, the 
growth intensity of grassland in the north of the Loess Plateau is not substantial, which further 
supports the expectation regarding to the protection of grassland and closed forest. In contrast, the 
increase of urban land pixels varies considerably under different scenarios, i.e., UES>HS>EPS, 
which possibly strengthens the conversion intensity in the surrounding areas and thus reduces 
crop production for the residents in these cities. 

The conversions of agricultural land, grassland, and urban land are noticeable (Fig. 4). Under 
all scenarios, urban land is predicted to expand and occupy agricultural land and other land use 
types. Grassland is predicted to increase in its area under both EPS and UES, and exhibits a 
trade-off relationship with deciduous broad-leaved forest but a synergic relationship with closed 
forest (Fig. S1). Additionally, closed forest shows the most rapid growth under EPS, and 
conversions between grassland and closed forest occur under all scenarios. Lastly, the others show 
trade-off relationships with grassland and urban land. The conversion of trade-off relationships 
among urban land, agricultural land, and ecological land (evergreen broad-leaved forest, 
deciduous broad-leaved forest, evergreen coniferous forest, closed forest, and grassland) to 
synergic relationships remains a complex challenge for regional land management. 
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Fig. 4 Conversions of land use types under HS (a), EPS (b), and UES (c) in the Loess Plateau during 2020-2050 


3.2 Variations in water retention services 


3.2.1 Spatial-temporal changes of precipitation and ETo 

In 2000-2015, both total precipitation and total ETo in the Loess Plateau changed smoothly (Fig. 5). 
The fluctuation range of precipitation was relatively large, and a negative correlation existed 
between the trend of precipitation and the trend of ETo (correlation coefficient of -0.673; P<0.05). 
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Fig. 5 Time series changes of pixel mean of precipitation and pixel mean of reference evapotranspiration (ETo) 
in the Loess Plateau from 2000 to 2015 


The average annual precipitation generally decreased from southeast to northwest in the Loess 
Plateau (Fig. 6a and b). The maximum values of 1197.15 mm in 2000 and 821.83 mm in 2015 were 
observed in the Wei-Fen River Basin, and the minimum values of 93.28 mm in 2000 and 95.56 mm 
in 2015 were found in the Lanzhou-Hekou zone. This variation was correlated with altitude. During 
2000-2015, the variation of precipitation was significant (Fig. 6c), and the areas with reduced 
precipitation extended from the middle to the east and west, especially in the Longyangxia-Lanzhou 
zone and Sanmenxia-Huayuankou zone, with the largest reduction of 412.77 mm. The areas with 
increased precipitation were mainly located in the north of the Loess Plateau, involving the 
Lanzhou-Hekou zone, Interior drainage zone, Hekou-Longmen zone, and north branch of the Haihe 
River Basin, and were also distributed in the south of the Wei-Fen River Basin, with a maximum 
increase of 170.67 mm. 

The average annual ETo was generally high in the north and low in the south of the Loess Plateau 
(Fig. 6d and e). The maximum values were 1200.89 mm in 2000 and 1150.72 mm in 2015, and the 
high values were concentrated in the Lanzhou-Hekou zone and Interior drainage zone. Obvious 
boundaries were observed between the areas with increasing and decreasing ETo changes (Fig. 6f), 
and the decrease mainly occurred in the Wei-Fen River Basin and south branch of the Haihe River 
Basin in the southeast of the Loess Plateau and the Lanzhou-Hekou zone in the northwest of the 
Loess Plateau, with a maximum decrease of 73.95 mm. 
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Fig. 6 Spatial distributions and changes of precipitation (a, b, c) and ETo (d, e, f) in the Loess Plateau from 2000 
to 2015 


3.2.2 Spatial-temporal changes of water retention services in 2000-2015 

In 2000-2015, the fluctuating changes of water yield and water retention were consistent in the 
Loess Plateau (Fig. 7), and the peak values occurred in 2003, 2007, and 2012, among which the 
peak value in 2003 was the highest (158.30x10° mm in the Loess Plateau). From 2000 to 2015, 
water yield decreased by 17.20x10° mm, while water retention increased by 0.09x10° mm, 
indicating that vegetation restoration measures such as the GGP exerted positive ecological 
effects. 
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Fig. 7 Time series changes of total water yield (a) and total water retention (b) in the Loess Plateau from 2000 
to 2015 


In 2000 and 2015, the spatial distribution of water yield was roughly the same, with the 
high-value areas mainly concentrated in the plains of the Wei-Fen River Basin, north branch of 
the Haihe River Basin, and south branch of the Haihe River Basin in the Loess Plateau (Fig. 8a 
and b). During this period (2000-2015), the area where water yield increased was mainly located 
in the north of the Loess Plateau (Fig. 8c), and the Interior drainage zone and its surrounding 
areas, which have experienced long-term large-scale vegetation restoration measures (Chen et al., 
2015) to convert arid deserts and deserts to oases in the recent years. However, water yield in the 
southern part mainly decreased. 

In addition, the spatial distribution of water retention has undergone major changes during 
2000-2015 (Fig. 8d—f), and there were significant changes in the north of the Interior drainage 
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zone and south of the Wei-Fen River Basin. Especially, high values of water retention were 
observed in the Interior drainage zone and its surrounding areas. Although the pixel mean of 
water retention in the Interior drainage zone increased by 1.93 mm, its value in the Loess Plateau 
decreased by 0.65 mm. 
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Fig. 8 Spatial distributions and changes of water yield (a, b, c) and water retention (d, e, f) in the Loess Plateau 
from 2000 to 2015 


3.2.3 Spatial-temporal changes of water retention services in 2020-2050 

In the simulation period (2020-2050), the multi-year average water yield in the Loess Plateau has 
the same spatial distribution patterns under two RCP scenarios (Fig. 9), and the high values are 
distributed in the southern plain of the Wei-Fen River Basin and the south branch of the Haihe 
River Basin. Water yield values are more prominent under RCP8.5 scenario, with a maximum 
pixel mean of 1050.00 mm. The spatial distributions of the multi-year average water retention 
values are different under two RCP scenarios (Fig. 10). Under RCP4.5 scenario, the south of the 
Wei-Fen River Basin, and the Interior drainage zone and its surroundings are the main areas with 
high values of water retention, with a maximum pixel mean of 277.46 mm. Under RCP8.5 
scenario, the high values are mainly located in the south of the Wei-Fen River Basin, with a 
maximum pixel mean of 330.20 mm. 

The pixel means of water yield under RCP4.5 and RCP8.5 scenarios are 96.63 and 95.46 mm, 
respectively. In the future (2020-2050), water yield values under RCP4.5 scenario comparatively 
resemble those in the historical period (2000-2015), and the distributions of areas with high and 
low water retention values under RCP4.5 scenario are similar to those in the historical period 
(Fig. 8). Under RCP8.5 scenario, the higher water retention values in the north of the Loess 
Plateau are considerably reduced and the total water retention decreases by 0.76x10° mm 
compared with that under RCP4.5 scenario. In addition, the percentage of pixel areas where water 
yield shows an increasing trend decreases with the order of RCP4.5-UES (71.57%)>RCP4.5-HS 
(71.33%)>RCP4.5-EPS (71.15%)>RCP8.5-UES (22.35%)>RCP8.5-HS (22.29%)>RCP8.5-EPS 
(22.22%); and the percentage of pixel areas where water retention shows an increasing trend 
decreases with the order of RCP4.5-HS (70.89%)>RCP4.5-EPS (70.88%)>RCP4.5-UES 
(70.78%)>RCP8.5-EPS (22.02%)>RCP8.5-HS (21.97%)>RCP8.5-UES (21.87%). For different 
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land use scenarios, there is a greater and more concentrated increase in water yield in the 
southeast of the Loess Plateau under UES (the standard deviations of the increased pixels for 
RCP4.5-UES and RCP8.5-UES are 80.97 and 170.28 mm, respectively), followed by HS (the 
standard deviations of the increased pixels for RCP4.5-HS and RCP8.5-HS are 61.46 and 143.97 
mm, respectively), and the increase in water yield under EPS is more balanced (the standard 
deviations of the increased pixels for RCP4.5-EPS and RCP8.5-EPS are 61.46 and 86.00 mm, 
respectively). 
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Fig. 9 Spatial distributions and changes of water yield under different scenarios (RCP4.5-HS, RCP4.5-EPS, 
RCP4.5-UES, RCP8.5-HS, RCP8.5-EPS, and RCP8.5-UES) in the Loess Plateau from 2015 to 2050 


At the plateau scale, total water yield values are as follows (in descending order): RCP4.5-UES 
(63.80x106 mm)>RCP8.5-UES (63.10x10° mm)>RCP4.5-HS (62.50x10° mm)>RCP8.5-HS 
(61.70x10° mm)>RCP4.5-EPS (61.90x10° mm)>RCP8.5-EPS (61.00x10° mm). Total water yield in 
the Loess Plateau is 0.76x10° mm higher under RCP4.5 scenario than under RCP8.5 scenario. The 
highest and lowest total water yield values are obtained under UES and EPS, respectively. The 
findings highlight the importance of increased vegetation cover in the conservation and fixation of 
water resources. The total water retention values at the plateau scale in descending order are as 
follows: RCP4.5-EPS (1.27x10° mm)>RCP4.5-HS (1.26x10° mm) and RCP4.5-UES (1.26x10° 
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mm)>RCP8.5-EPS (0.91x10° mm)>RCP8.5-HS (0.90x10° mm)>RCP8.5-UES (0.89x10° mm). 
Total water retention is 0.37x10° mm higher (i.e., more favourable for water retention) under 
RCP4.5 than under RCP8.5. 
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Fig. 10 Spatial distributions and changes of water retention under different scenarios (RCP4.5-HS, RCP4.5-EPS, 
RCP4.5-UES, RCP8.5-HS, RCP8.5-EPS, and RCP8.5-UES) in the Loess Plateau from 2015 to 2050 


3.2.4 Importance rankings of water retention services 


The importance of water retention services in the Loess Plateau in the historical period 
(2000-2015) and in the future (2020-2050) under RCP4.5 and RCP8.5 scenarios can be divided 
into five levels (general (0.00-3.00 mm), more (3.00—10.00 mm), moderate (10.00-30.00 mm), 
high (30.00-80.00 mm), and extreme (80.00—330.20 mm)) using natural breaks method, and the 
percentages of areas with these different levels were calculated (Table 1). Areas with moderate to 
high importance levels are mainly located in the northwest of the Lanzhou-Hekou zone, Interior 
drainage zone, Hekou-Longmen zone, and southern Wei-Fen River Basin (Fig. 11). Under 
RCP4.5 scenario, 4.33% functional areas have moderate to high importance levels of water 
retention services, which is 2.23% greater than that under RCP8.5 scenario. Further, under 
RCP8.5 scenario, the importance levels of water retention services show decreasing trend in the 
northwest of the Loess Plateau and increasing trend in the south of the Loess Plateau. Functional 
areas with moderate, high, and extreme importance levels only account for 2.10% of the total area 
under RCP8.5 scenario. These ranking results may provide reference data for the exploitation of 
water resources and the determination of priority points and key areas for the enhancement of 
water retention in the Loess Plateau. 
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Table 1 Percentage of area with different importance levels of water retention services in the Loess Plateau in 
the historical period (2000-2015) and in the future (2020—2050) under RCP4.5 and RCP8.5 scenarios 


Percentage of area with different importance levels of water retention services (%) 


Period 


General More Moderate High Extreme 
(0-3 mm) (3-10 mm) (10-30 mm) (30-80 mm) (80—330.2 mm) 
Historical period (2000-2015) 81.29 12.92 5.20 0.59 0.05 
Future (2020-2050) under 
RCP45 84.65 11.02 4.03 0.29 0.01 
Future (2020-2050) under 
RCP85 88.54 9.36 1.91 0.18 0.01 


Note: RCP, representative concentration pathway. The importance of water retention services in the Loess Plateau in the historical 
period (2000-2015) and in the future (2020-2050) under RCP4.5 and RCP8.5 scenarios can be divided into five levels using natural 
breaks method: general, 0.00-3.00 mm; more, 3.00-10.00 mm; moderate, 10.00-30.00 mm; high, 30.00-80.00 mm; and extreme, 
80.00—330.20 mm. 
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Fig. 11 Spatial distributions of different importance levels of water retention services in the Loess Plateau in the 
historical period (2000-2015; a) and in the future (2020-2050; b and c) under RCP4.5 and RCP8.5 scenarios. The 
importance of water retention services in the Loess Plateau in the historical period (2000-2015) and in the future 
(2020-2050) under RCP4.5 and RCP8.5 scenarios can be divided into five levels using natural breaks method: 
general, 0.00-3.00 mm; more, 3.00—10.00 mm; moderate, 10.00-30.00 mm; high, 30.00-80.00 mm; and extreme, 
80.00—330.20 mm. 


4 Discussion 


4.1 Correlation between water retention and water yield 


Although large-scale vegetation restoration has been carried out in the Loess Plateau, the total 
amount of water retention is still low in the whole region, and rain retention capacity is very 
finite. Water retention in the Loess Plateau only accounts for 2.11% of water yield, while water 
yield accounts for about 24.25% of precipitation. After vegetation restoration, increased water 
consumption and canopy water consumption consume high amounts of water resources in the 
Loess Plateau. Trying to obtain high water yield and low water consumption through vegetation 
restoration and other measures is a challenge. Sun et al. (2006) showed that afforestation 
measures under the average historical climate conditions in China would decrease runoff at rates 
of 50.00—-150.00 mm/a or with percentages of 20.00%—40.00% in most areas of the country, 
especially in the Loess Plateau (decreasing rate of 50.00 mm/a or decreasing percentage of 
50.00%). According to Lü et al. (2012), after the implementation of the GGP (2000-2008) in the 
Loess Plateau, runoff in the region decreased at a rate of 10.30 mm/a, and only less than 10.00% 
of the whole region had increased water yield. In addition, according to Feng et al. (2016a), the 
massive increase in vegetation after the implementation of the GGP increased net primary 
productivity (NPP) and ET» in the Loess Plateau, and further vegetation expansion could cause 
regional water shortage. Furthermore, Shao et al. (2019) reported that large-scale vegetation 
restoration under the implementation of the GGP increased ETo significantly (4.39 mm/a) in the 
Loess Plateau. In 2001-2015, the increase of vegetation under the implementation of the GGP 
resulted in the consumption of approximately 31.00x10* m?/a ETo. Wang et al. (2021) pointed out 
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that water retention in the Zhanghe River Basin in the eastern part of the Loess Plateau decreased 
by 1.10x10® m*/decade from 1960 to 2016, mainly due to increased ecological engineering 
measures and reduced water consumption by human activities. 

Consistent with the results of Feng et al. (2016a), in this study, ETo was significantly negatively 
correlated with precipitation, water yield, and water retention, with correlation coefficients of 
—0.673, —0.720, and —0.737, respectively (P<0.01). In the Loess Plateau, the areas with high rates 
of water retention to water yield were mainly located in the Lanzhou-Hekou zone, Interior 
drainage zone, and Hekou-Longmen zone in the northwest, with values of 2.58%, 6.25%, and 
1.96%, respectively. The complex influences of ecological projects (e.g., the GGP) and the spatial 
heterogeneity of regions should be comprehensively considered to enable the formulation of 
suitable ecological measures in the Loess Plateau. 


4.2 Variations in water retention services in a typical sandy land area in the Loess Plateau 


Increased vegetation coverage in the Loess Plateau has effectively controlled the sediment 
transport into the Yellow River (Fu et al., 2017). The Mu Us Sandy Land is a typical example; 
it is one of the major sandy lands in China and is located in the north of central Loess Plateau 
(Fig. 12). Its harsh ecological environment limits socio-economic development; however, after 
the implementation of various ecological and environmental governance measures for over 70 
a, more than 90.00% of the desertified land is under management (Li et al., 2017; Liu et al., 
2020). From 2000 to 2015, 91.75% of area in the Mu Us Sandy Land has increased in water 
retention, and the total amount of water retention has significantly increased from 0.16x10° to 
0.37x10° mm. 
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(c) Water retention under RCR8.5-EPS in 2050 
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Fig. 12 Spatial distributions of water retention change in the Mu Us Sandy Land in 2000-2015 (a), and spatial 
distributions of predicted water retention in the Mu Us Sandy Land in 2050 under RCP4.5-EPS (b) and 
RCP8.5-EPS (c) 
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Vegetation restoration projects are the key factors that enhance regional water retention 
services (Liu et al., 2020) and can significantly increase the ecological, social, and economic 
benefits to the Loess Plateau. Under RCP4.5-EPS, vegetation in the Mu Us Sandy Land is 
expected to increase continuously, thus considerably promoting water retention. In addition, the 
total water retention in the Mu Us Sandy Land in 2050 is expected to reach 0.46x10° mm, which 
is much higher than that under RCP8.5-EPS (0.13x10° mm). The series of long-term ecological 
restoration projects have significantly improved the typical ecological environment and enhanced 
water retention capacity of the Loess Plateau. 


4.3 Uncertainty analysis 


Integrated modelling methods can facilitate the formulation of land use planning policies that 
consider both socio-economic development and ecological conservation activities. In this study, 
we conducted the multi-scenario simulation of water retention services by integrating the LUCC 
models, water retention assessment, and the spatial mapping methods of intermediate parameters 
(such as Zhang coefficient), split the process "black box" of mainstream models such as InVEST, 
and improved the limitation of its constant settings to better depict the actual conditions of the 
Loess Plateau across multiple climatic zones. Specifically, the improved Markov model was used 
to predict the land use demand linearly, but in reality, the threshold and inflection point of 
non-linear change of land use should be taken into consideration, and the simulation based on the 
FLUS model should further consider the role of policy and so on. Although water retention 
service assessment was processed to provide an interface for the verification of process 
parameters, the lack of verification datasets needs to be resolved urgently. 

Simulations under six scenarios demonstrated that water yield and water retention in the Loess 
Plateau are highly consistent with precipitation fluctuations; the average correlation coefficients 
between water yield and precipitation and between water retention and precipitation are 0.93 and 
0.91, respectively, which are consistent with the results of previous studies (Feng et al., 2012; Su 
and Fu, 2013; Fu et al., 2017). The results of simulation indicate that climate change has a direct 
and significant impact on water retention services in the Loess Plateau. Moreover, the impact of 
LUCC and trade-offs cannot be ignored, especially the LUCC caused by two major human 
activities, i.e., urbanization and vegetation restoration. It is to be noted that, in this study, the 
forecast of land use demand under different land use scenarios employs a linear calculation 
process. It can be used as an effort on revealing the greatest impact of different land use/land 
cover patterns through this more "extreme" forecast, practically along two main development 
routes of current and future ecological protection and urban land expansion. In addition, 
quantifying the impacts of LUCC and climate change is an important extension of this study. 


5 Conclusions 


In the Loess Plateau, the main characteristics of LUCC in the historical period (2000-2015) and 
in the future (2020-2050) are the mutual conversions among agricultural land, urban land, and 
grassland. There is a significant trade-off between the expansion of urban land and the decrease of 
agricultural land. It is thus recommended to strictly control the expansion of cities and towns and 
to strengthen the conservation and quality improvement of agricultural land and ecological land. 
In the historical period (2000-2015), the change trends of water yield and water retention are 
consistent in the Loess Plateau, and the increase of water retention in the Interior drainage zone 
reflects the effectiveness of vegetation restoration. In the future (2020-2050), the simulated 
spatial patterns of water yield under two RCP scenarios (RCP4.5 and RCP8.5) are similar in the 
Loess Plateau, while RCP4.5 model is more conducive to water retention. Compared with HS and 
UES, water retention under EPS is relatively high. Based on the above findings, we propose that 
the Loess Plateau should actively explore a new development model that focuses on ecological 
protection with equal emphasis on quality and benefits on the basis of land intensive 
development. 
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Appendix 


Table S1 Data sources and data processing used in this study 


Name Source Format Period Processing method 
LULC European Space Agency ; 300 m, 2000-2015 Reclassification and 
(https://www.esa-landcover-cci.org) grid resampling 
Resource and Environment Data Cloud Platform, the 1000 m Slope, aspect, and 
DEM Chinese Academy of Sciences id ? 2003 catchment calculations 
(http://www.resdc.cn/Default.aspx) gm by ArcGIS 10.2 
Soil properties Environmental and Ecological Science Data Center for 30s, 2013 Projection conversion 
West China (http://westdc.westgis.ac.cn/) grid and resampling 
China Meteorological Data Service Center Text 2000-2015 Data organization and 
Precipitation (http://data.cma.cn/) E , interpolation , 
National Aeronautics and Space Administration (NASA) 25 km, 2020-2050 Format conversion and 
(https://nex.nasa.gov/nex/projects/1356/) grid resampling 
Meteorological China Meteorological Data Service Center Data organization and 
elements (http://data.cma.cn/) Tert 200072019 interpolation 
Land-Atmosphere Interaction Research Group at Sun 30s Projection and 
LAI Yat-sen University, China a 2000-2015 ede in 
(http://globalchange.bnu.edu.cn/) gn pung 
Resource and Environment Data Cloud Platform, the 1000 m 
Population, GDP Chinese Academy of Sciences grid a 2000-2015 Projection 


(http://www.resdc.cn/Default.aspx) 
Resource and Environment Data Cloud Platform, the 


ae Chinese Academy of Sciences Vector laa andara 
(http://www.resdc.cn/Default.aspx) pping 
Catchment Resource and Environmental Science Data Center, the Projection and data 
: Chinese Academy of Sciences Vector f 
boundaries cropping 


(http://www.resdc.cn/data.aspx) 
Note: LULC, land use/land cover; DEM, digital elevation model; LAI, leaf area index; GDP, gross domestic product. 
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Fig. S1 Correlation of different land use types under HS (a), EPS (b), and UES (c) in the Loess Plateau. HS, 
historical scenario; EPS, ecological protection scenario; UES, urban expansion scenario. 


